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Abstract 

Using nearly 60 pb _1 of data collected with the CLEO-c detector at the ^(3770) resonance, 
we measure absolute branching fractions for three D° and two D + Cabibbo-allowed hadronic 
decay modes, and the cross section for e + e~ — > DD at E cm = 3.77 GeV. We report preliminary 
measurements of the reference branching fractions, B(D° — > R~ir + ) = (3.92 ± 0.08 ± 0.23)% and 
B(D + — > K~it + it + ) = (9.8 ±0.4 ±0.8)%, and preliminary measurements of other major branching 
fractions, B(D° -> R-tt+tt ) = (14.3 ± 0.3 ± 1.0)%, B(D° -> R-it+it+it-) = (8.1 ± 0.2 ± 0.9)%, 
and B(D + — > ivT^7r + ) = (1.61 ± 0.08 ± 0.15)%. We determine preliminary values of the cross 
sections, a(e+e~ -» D°Z) ) = (3.47 ± 0.07± 0.15) nb, a(e+e- -» D+D~) = (2.59 ± 0.11 ± 0.11) nb, 
and a(e + e~ — > Z)Z)) = (6.06 ± 0.13 ± 0.22) nb. We note that the Monte Carlo simulations used 
in calculating efficiencies in this analysis included final state radiation. However, the branching 
fractions used in the Particle Data Group averages do not include this effect. If we had not included 
final state radiation in our simulations, branching fractions would have been 0.5% to 2% lower. 



*Submitted to the 32 nd International Conference on High Energy Physics, Aug 2004, Beijing 
'''On leave of absence from University of Chicago. 



2 



I. INTRODUCTION 



We present preliminary absolute measurements of the Cabibbo-allowed D° and D + 
branching fractions, B(D° -> K~ti + ), B(D° -> K'n + n°), B(D° -> ^-tt+tt+tt-), £(£>+ -> 
A'~7r + 7r + ), and B(D + — ► i^°7r + ). Two of these branching fractions, B(D° — > K~tt + ) and 
B(D + — ► i^ _ 7T + 7r + ), along with the branching fraction for — ► 07r + , are particularly 
important because essentially all other £)°, D + , and Z)+ branching fractions are determined 
from ratios to one or the other of these branching fractions As a result, nearly all 
branching fractions in the weak decay of heavy quarks are ultimately tied to one of these 
three branching fractions, called reference branching fractions in this paper. Furthermore, 
these reference branching fractions appear in many measurements of CKM matrix elements 
for c and b quark decay. 

We note that the Monte Carlo simulations used in calculating efficiencies in this analysis 
included final state radiation. However, the branching fractions used in the Particle Data 
Group averages do not include this effect. If we had not included final state radiation in our 
simulations, the branching fractions we measure would have been 0.5% to 2% lower than 
they are. 



II. BRANCHING FRACTIONS AND PRODUCTION CROSS SECTIONS 

The data for these measurements were obtained in e + e~ collisions at E cm = 3.77 GeV, 
the peak of the ^(3770) resonance. At this energy, no additional hadrons accompany the 
D°D° and D + D~ pairs that are produced. These unique DD final states provide a powerful 
tool for addressing the most vexing problem in measuring absolute D branching fractions at 
higher energies - the difficulty of accurately determining the number of D mesons produced. 
Following the MARK III collaboration 0, 01 , reconstruction of one D meson (called single 
tag or ST) serves to tag the event as either D°D° or D + D~. Branching fractions for 
D° or D + decay can then be obtained from the number of "double tag" (DT) events in 
which both the D and the D are reconstructed, without knowledge of luminosity or the 
total number of DD events produced. Also, since E cm is below threshold for production of 
DfD~, accurate absolute measurement of B(D^ — > 07r + ) using the techniques in this paper 
remains a challenge for the future. 

If no distinctions are made between D and D decays and their efficiencies, the number 
of ST events Ni observed in the decay mode % (and its charge conjugate) with branching 
fraction £>j and detection efficiency e, will be, 

Ni = 2N DD Bie h 

where N D jj is the number of DD events produced in the experiment. Then, the number of 
DT events with DD pairs reconstructed in modes % and j will be, 

N ij = 2N DD B i B J t ij (ifz^j) or N ii = N DD B 2 i e ii (if i = j) 

Hence, the ratios of DT events (Nj) to ST events (Nj) provide absolute measurements of 
the branching fraction Ej, 

Bi = ~T7~ or Bi = — . 

Nj eij N en 
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Note that ~ SiSj, so ej nearly cancels from the ratio ej/eij. Hence, a branching fraction 
Bi obtained using this procedure is nearly independent of the efficiencies of tagging modes. 
Of course, Bi is sensitive to e« and its uncertainties. 

Estimating errors and combining measurements using just these expressions is very diffi- 
cult because iV^ and Nj are correlated (whether or not i = j) and measurements of Bi using 
different tagging modes j are also correlated. To address this problem, we have developed a 
X 2 fitting procedure Q which fits simultaneously all charged and neutral D branching frac- 
tions and the numbers of charged and neutral DD pairs. Although the D° and D + branching 
fractions are statistically independent, systematic effects introduce significant correlations 
among them. Therefore, we fit both charged and neutral D mesons simultaneously, and 
we include both statistical and systematic uncertainties, as well as their correlations, in the 
fit. We also perform corrections for backgrounds, efficiency, and crossfeed among modes 
directly in the fit, as the sizes of these adjustments depend on the fit parameters. Thus, 
all experimental measurements, such as event yields, efficiencies, and background branching 
fractions, are treated in a consistent manner. We actually fit D and D yields separately 
in order to accommodate possible differences in efficiency, but charge conjugate branching 
fractions are constrained to be equal. 

The number of DD events that were produced can be calculated from 

N - - MiliL or N - — 

The production cross sections for D°D° and D + D~ can then be obtained by combining 
N D 0£)0 and Nu+ £>- - which are determined in the branching fraction fit - with the integrated 
luminosity C. Note that N D p obtained by this procedure is almost independent of the 
efficiencies. 

Charge conjugate particles and decay modes are always implied in this paper unless stated 
otherwise. 

III. THE CLEO-C DETECTOR 

The CLEO-c detector is a modification of the CLEO III detector 0, S Q m which the 
silicon-strip vertex detector was replaced with a six-layer vertex drift chamber, whose wires 
are all at small stereo angles to the beam axis The charged particle tracking system, 
consisting of the vertex drift chamber and a 47-layer central drift chamber [f| operates in 
a 1.0 T magnetic field, whose direction is along the beam axis. The momentum resolution 
achieved with the tracking system is ~ 0.6% at p = 1 GeV/c. Photons are detected in an 
electromagnetic calorimeter consisting of about 8000 CsI(Tl) crystals [Hf. The calorimeter 
attains a photon energy resolution of 2.2% at E 1 — 1 GeV and 5% at 100 MeV. The solid 
angle coverage for charged and neutral particles of the CLEO-c detector is 93% of 4ir. We 
utilize two devices to obtain particle identification information to separate K ± from tt ± : 
the central drift chamber (which provides measurements of ionization energy loss - dE/dx) 
and a cylindrical ring-imaging Cherenkov (RICH) detector surrounding the central drift 
chamber. The solid angle of the RICH detector for separation of K ± from 7r =1= is 80% of 47r. 
The combined dE/dx-TtiCH particle identification procedure has a pion or kaon efficiency 
> 90% and a probability of pions faking kaons (or vice versa) < 5%. 

The integrated luminosity (£) of the data sets was measured using 77 events in the 
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calorimeter [9|. Event counts were normalized with a Monte Carlo (MC) simulation based 
on the Babayaga event generator combined with GEANT-based detector modeling. 

IV. DATA SAMPLE AND EVENT SELECTION 

In this analysis, we utilized a total integrated luminosity of C = 57.2 ± 1.7 pb _1 of 
e + e _ data collected at E cm = 3.77 GeV. The data were produced with the Cornell Electron 
Storage Ring (CESR) operating in a new configuration |8j that includes six wiggler magnets 
to enhance synchrotron radiation damping at energies in the charm threshold region. Since 
these data were obtained, six additional wiggler magnets have been installed in CESR and 
future CESR operation in this energy region will be with the full complement of twelve 
wiggler magnets. The spread in E cm with the six wiggler magnets is 2.3 MeV. 

We reconstructed all events that satisfied the CLEO-c trigger and preserved them for 
further analyses. Most CLEO-c analyses of D meson decays utilize a subset of the data 
in which each event contains at least one D meson candidate selected by a standard set of 
requirements. This selection procedure begins with standardized requirements for it^, K ± , 
7T°, and Kg candidates. 

Charged track candidates must pass the following quality requirements: the momentum 
of the track p must be in the range 0.050 < p < 2.0 GeV/c, the track must pass within 
0.5 cm of the origin in the x-y plane (transverse to the beam direction), the track must 
pass within 5.0 cm of the origin in the z direction, the polar angle 9 was required to be in 
the range | cos#| < 0.93, and the number of hits reconstructed on the track was required to 
be at least half of the number of layers traversed by the track. Both requirements on the 
distance of the track from the origin of the coordinate system (which is close to the point 
where the two beams intersect) are approximately five times the standard deviation for the 
corresponding measurement. 

We identified charged track candidates as pions or kaons using dE/dx and RICH informa- 
tion. Each track was considered as both a potential K and 7r ± candidate. In the rare case 
that no useful information of either sort was available, we did not change this assignment. 

If dE/dx information was available, we calculate an d Xe(K) horn the dE/dx mea- 

surements, the expected dE/dx for pions and kaons of that momentum, and the appropriate 
standard deviations. We rejected tracks as kaon candidates when Xe{K) was greater than 
9, and similarly for pions. The difference xM 71 ") — Xe{K) was a l so calculated. If dE/dx 
information was not available this x 2 difference was set equal to 0. 

We used RICH information if the track momentum was sufficiently above the K threshold 
(p > 0.55 GeV/c), and the track was within the RICH acceptance (| cos^l < 0.8). We then 
rejected tracks as kaon candidates when the number of Cherenkov photons detected for the 
the kaon hypothesis was less than three, and similarly for pions. When there were at least 
three photons for the hypothesis under consideration, and information was available for both 
pion and kaon hypotheses, we calculated a x 2 difference for the RICH, x'lti 71 ') ~Xr{K)i from 
the locations of Cherenkov photons and the track parameters. Otherwise, we set this x 2 
difference equal to 0. 

Then we combined these x 2 differences in an overall x 2 difference, A% 2 = x%(^)~ Xe(K) + 
Xr^) — Xr{K). If A% 2 < 0, we used the track as a pion candidate; If Ax 2 > 0, we used 
the track kaon candidate. 

We formed neutral pion candidates from pairs of photons, each of whose energy was 
greater than 30 MeV and whose showers pass photon quality requirements. An uncon- 
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strained mass M (77) was calculated from the energies and momenta of the two photons. 
This mass was required to be within 3a with a nominal 7r° mass value which varied slightly 
with the total momentum of the 7r° candidate. The energy and momentum obtained from 
a kinematic fit of the two photon candidates to the mass M n o from the PDG was then 
used in further analysis. 

We built Kg candidates from pairs of tracks that intersect a vertex. We then subjected 
the tracks to a constrained vertex fit and used the resulting track parameters to calculate 
the invariant mass, M(7r + 7r~). We called the track pair a K s candidate if the invariant mass 
M(tt + ti-) was within 3a, with a = 4 MeV/c 2 , of the mass M K o from the PDG Q. 

After event reconstruction, the next major step was selecting DD events with relatively 
loose requirements for further analysis. In this CLEO-c standard selection process we re- 
constructed D° candidates in 18 important decay modes and D + candidates in 6 important 
decay modes, although only three D° modes and two D + modes were used in this analy- 
sis. Two requirements were placed on D° and D + candidates formed from 7r ± , K ± , 7r° and 
K s candidates selected using the requirements described above. We calculate the energy 
difference, AE = E(D) — E , where E(D) is the total energy of the particles in the D 
candidate and Eq is the energy of the e beams. We accept candidates for further analysis 
if |Ai?| < 100 MeV. We calculate the mass of the D candidate M(D) by substituting the 
beam energy E for the energy E(D) of the D candidate, i.e., M 2 (D)c 4 = E 2 — p 2 (D)c 2 , 
where p(D) is the total momentum of the particles in the D candidate. Then we required 
M(D) > 1.83 GeV/c 2 to retain a D candidate for further analysis. We included all events 
with D candidates in any of the 24 decay modes that satisfied these loose requirements. 

For the analysis described in this paper, we restricted our attention to the five decay 
modes mentioned previously, and we refined the requirements for acceptable D candidates 
that were described in the previous paragraph. First, we used the more restrictive require- 
ments on IA.EI, given in Table [Q that were tailored for each individual decay mode. For 
the ST analysis, we chose the candidate with the smallest |Ai?|, if there was more than 
one D candidate in a particular mode. Multiple candidates were rare in all modes except 
D° — > K~7t + 7t where approximately 15% of the events had multiple candidates. In events 
with only two charged tracks that were consistent with our requirements for D° —>■ K~n + de- 
cays, we also imposed additional requirements to eliminate e + e~ — > e + e~ and e + e~ — > fi + fi~ 
events. Since we are not considering any all-neutral modes in this analysis, these require- 
ments only affect ST yields, i.e., we would not try to find DT candidates with the D or D 
decaying in an all neutral mode recoiling from a D° — > K~tt + candidate. 

TABLE I: Requirements on AE for D candidates accepted for the final analysis. The limits are 
set at approximately 3a. 



Mode 




Requirement (MeV) 


D° — > 




\AE\ < 29.4 


D° — > 


K-ir + ir° 


-58.3 < AE < 35.0 


D° — > 




\AE\ < 20.0 


D+ -. 




|A£7| < 21.8 


D+ — : 


■KgTT+ 


\AE\ < 26.5 
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To obtain DD events for DT yields, we select only one candidate per event per combina- 
tion D and D decay modes. We apply the AE requirements described for the ST analyses, 
but do not choose the best candidate on the basis of minimum |A£J|. Instead, we choose 
the combination with the average of M(D) and M(D) - i.e., M(D) = [M (D) + M(D)]/2 - 
closest to M D . In careful studies of Monte Carlo events, we demonstrated that this proce- 
dure does not generate false peaks at the D mass in the M(D) vs. M(D) distributions that 
are narrow enough or large enough to be confused with the DT signal. 



V. GENERATION AND STUDY OF MONTE CARLO EVENTS 

We used Monte Carlo simulations to develop the procedures for measuring branching 
fractions and production cross sections, to understand the response of the CLEO-c detector, 
to determine parameters to use in fits to determine data yields, and to estimate and under- 
stand possible backgrounds. In each case e + e~ — > ?/>(3770) — > DD events were generated 
with the EvtGen program 12 . and the response of the detector to the daughters of the DD 



decays was simulated with GEANT The EvtGen program includes simulation of initial 
state radiation (ISR) - radiation of a photon by the e + and/or e~ before their annihilation. 
The program PHOTOS 0] was used to simulate final state radiation (FSR) - radiation 
of photons by the charged particles in the final state. We simulated two types of D meson 
decays: 

• signal Monte Carlo decays, in which a D or a D decays in one of the 5 modes measured 
in this analysis, and 

• generic Monte Carlo decays, in which a D or D decays in accord with d ecay modes 
and branching fractions based on the 2002 Particle Data Group compilation [1J| . (Some 
tuning was required to match inclusive distributions in data.) 

Using these two types of simulated decays, we generated three types of Monte Carlo events: 

• generic Monte Carlo events, in which both the D and the D decay generically, 

• single tag signal Monte Carlo events, in which either the D or the D always decays in 
one of the 5 modes measured in this analysis while the D or D, respectively, decays 
generically, and 

• double tag signal Monte Carlo events, in which both the D and the D decay in one of 
the 5 modes in this analysis. 

Our signal Monte Carlo events included simulations of ISR, but our generic Monte Carlo 
events were generated before ISR was included in the simulations. The roles played by these 
Monte Carlo samples are described below. 

We applied the same selection criteria for D candidates and DD events in analyzing 
data and Monte Carlo events. We compared many distributions of particle kinematical 
quantities in data and Monte Carlo events to asses the accuracy and reliability of the Monte 
Carlo simulation of event generation and detector response. The agreement between data 
and Monte Carlo events for both charged and neutral particles was excellent for nearly all 
distributions of kinematic variables that we studied. The results of this analysis are not 
sensitive to the modest discrepancies that were observed in a few distributions. 
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VI. SINGLE TAG EFFICIENCIES AND DATA YIELDS 



We used binned likelihood fits to ST M(D) distributions in Monte Carlo events to deter- 
mine experimental resolutions and ST efficiencies. We then used some of these experimental 
resolution parameters in binned likelihood fits to data to determine ST data yields. Four 
probability distribution functions were used in fits to ST data and ST Monte Carlo events: 

• Fits to mass (m = M(D)) distributions in ST Monte Carlo events without ISR used a 
relatively narrow core Gaussian g(m; pi, cr) to account for beam energy spread and a small 
contribution from detector resolution. The mean p, standard deviation cr, and area A of 
the core Gaussian were among the parameters determined in the fits. 

• In fits to Monte Carlo events with ISR and to data, the core Gaussian function was 
replaced by an inverted Crystal Ball |l6j function c(m; p, a, a, n) to account for the ra- 
diative tail on the high mass side of the peak due to ISR. This function is identical to 
g(m; p, cr) for m < p + aa, but for larger values of m the Gaussian is replaced with a 
radiative tail function that depends on a and n. The parameters a and n are determined 
by fits to the mass distributions. 

• All fits included a wider bifurcated Gaussian, b(m; p, ol, <tr) with different standard 
deviations <jl and <tr to the left and right of the peak at m = p, respectively. In all 
fits, we constrained the parameter p in this function to be equal to the mean of the core 
Gaussian. This term models misreconstruction of charged particle tracks, neutral pions, 
and neutral kaons. The ratios p^ = ctl/ct, p R = o-r/ct of the left and right standard 
deviations of the bifurcated Gaussian to the standard deviation of the core Gaussian and 
the ratio (pa) of the area of the bifurcated Gaussian to the total signal area (i.e., area 
of the bifurcated Gaussian plus the area of either the core Gaussian or the Crystal Ball 
function) were parameters in the fits. In fits involving ISR, the parameters p and a used 
here refer to the Gaussian parameters of the Crystal Ball function. 

• Although combinatorial backgrounds in all M(D) distributions were very small, an AR- 
GUS background function [l5j included in each of the fits to account for this contribution. 
The ARGUS function a(m; Eq, £) is a function of the beam energy Eq and a parameter £ 
determined in the fit. 

We always fit the M(D) and M(D) distributions for ST events in each signal mode indi- 
vidually, with the parameters of the signal distributions constrained to be equal, but with 
no constraints on background sizes or shapes. The parameters determined from M(D) and 
M(D) distributions for a given mode were in excellent agreement. The manners in which 
we utilized these four contributions to fit Monte Carlo events and data are: 

• Single tag yields in signal Monte Carlo events without ISR 

We obtained signal Monte Carlo events without ISR from our sample of signal Monte 
Carlo events with ISR by rejecting events with ISR photons with energies E 1 > 25 keV. 
We used this sample to study the shapes of the experimental M(D) distributions due 
to detector resolution and misreconstructions and determine the parameters of the core 
and bifurcated Gaussians for use in fits to other distributions. The means of the core 
Gaussians are all in excellent agreement with the masses M^o and Mu+ from the PDG £0] . 
The widths of the core Gaussians are in the range a = 1.25-1.35 MeV for all modes except 
D° — > K~7t + 7t° where cr ~ 1.7 MeV. The contributions of the bifurcated Gaussians are 
small and their widths are somewhat larger than the width of the core Gaussian: for all 
modes except D° — > A"~7r + 7r° we find pi = 2.0, Pr = 2.5-2.6, and pa = 3%-6%; while 
for D° — » A"~7T + 7r° we find pi = 2.9, pr = 3.2, and pa = 11%. The differences between 
D° — > K~7T + 7T and the other modes are due to 7r° reconstruction uncertainties. 
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Single tag yields in signal Monte Carlo events with ISR 

We used all single tag Monte Carlo events to determine the efficiencies £j for detecting 
ST decays in the signal modes and to determine the contribution of ISR to the shape of 
the M(D) peak. For each signal mode, we used the Gaussian parameters, pl, and 
Pa determined from the ST events without ISR. We used the total areas (sum of the two 
Gaussians and the radiative tail) of the M(D) peaks to determine the ST efficiencies e^. 
Single tag yields in generic Monte Carlo 

We used our large sample of generic Monte Carlo data (without ISR) to verify the ef- 
fectiveness of our fitting procedures, to study combinatorial backgrounds, and to provide 
input data to test the branching fraction fit program. We fixed the values of the pa- 
rameters a, Pl, and pR to the values obtained from signal Monte Carlo events without 
ISR. 

Single tag yields in data 

We fit the M(D) distribution for a given decay mode in ST data using the Gaussian 
parameters, pi, Pr, and p\ determined from the signal Monte Carlo events without ISR. 
From the fit, we determined a, A, a, n, and £. The yields are the sums of the contributions 
from the two Gaussians and the radiative tail. The data and the results of these fits are 
illustrated in Figs. and El for D°D° and D + D~ candidate events, respectively. The 
linear plots demonstrate that the misreconstruction contributions and the combinatorial 
backgrounds are small, and the semi-logarithmic plots demonstrate that the signals and 
backgrounds are modeled well by the fits. 

TABLE II: Single tag data yields and efficiencies and their statistical uncertainties. 
D or D Mode Yield (10 3 ) Efficiency 



D° — > 






5.14 ±0. 


.07 


65.1 ±0.6 


D° —> 


K + it- 




5.16 ±0. 


.08 


66.3 ±0.6 


D° — > 






9.62 ± 0. 


.12 


33.6 ±0.4 


D° — > 


K+TT- 


-7T° 


9.58 ±0. 


.12 


34.0 ±0.4 


D° — > 




- + — 
7T 7T 


7.39 ± 0. 


.10 


45.1 ±0.5 


D° — > 


K+TT~ 




7.39 ± 0. 


.10 


45.5 ±0.5 


D+ -> 


■ K--K- 




7.58 ±0. 


09 


52.2 ±0.5 


D~ -. 


■ K+TT- 




7.57 ±0. 


09 


51.9 ±0.5 


D+ -> 






1.09 ±0. 


.04 


45.6 ±0.5 


D~ — 3 


■K° S 7T- 




1.12 ±0. 


.04 


45.9 ±0.5 



Table |H] gives the ST yields in data and the ST efficiencies calculated from the signal 
Monte Carlo events with ISR. These quantities are used in the fit, described in Sec. El for 
branching fractions and numbers of DD events. 

VII. DOUBLE TAG EFFICIENCIES AND DATA YIELDS 

The components that must be included in the DT fits are illustrated in Fig. El which 
shows distributions of M(D) vs. M(D) for DT event candidates that are reconstructed in 
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the D° — > K + 7t~7t° and _D° — > i^~7r + 7r° modes. Figure 0(a) shows all candidates (including 
multiple candidates) after imposition of the AE requirement described in Sec. IIVI The 
principal features of this two-dimensional distribution are: 

• There is an obvious strong signal concentration with a complicated shape in the region 
surrounding M(D) = M{D) = M D o. 

• There is a radiative tail above the signal peak at 45° to the axes. This correlation is 
due to the fact that - neglecting measurement and reconstruction errors - the values of 
M(D) and M(D) calculated using the beam energy will both to be too large by the same 
amount if energy was lost due to ISR. 

• There are faint horizontal and vertical bands at M(D) = M D o and M(D) = M D o, 
respectively. The events in these bands are events in which the D° {D°) candidate was 
reconstructed correctly, but the D° (D°) was misreconstructed. 

• There is a somewhat more visible band below the peak at 45° to the axes that continues 
through the signal region and the radiative tail. This band is populated by two sources 
of background: 

o There are DD events, in which all of the particles were found and reconstructed 
reasonably accurately, but one or more particles from the real D° were interchanged 
with the corresponding particles from the D°, e.g., the two 7r°s were interchanged. 

o There are also continuum events, i.e., annihilations into uu, dd, and ss quark pairs. 
As illustrated in Fig. 0(b), all three bands below the peak are depopulated noticeably after 
choosing the candidate with the average mass M(D) nearest to M^o. Since careful studies 
of this choice using Monte Carlo simulations of generic DD decays and continuum events 
showed no evidence of false peaks in the region of the D mass, the net effect of this procedure 
is reduction of background without introduction of a false signal. 

We included five terms in unbinned likelihood fits to account for these features of the 
two-dimensional M(D) vs. M(D) distributions. These terms were functions of m — M(D), 
m = M(D), m = M(D) = [M(D) + M(D)}/2, and 5m = [M(D) - M(D)}/2. The terms 
are: 

• c(m; p, a, a, n) g(5m; p,, as), a Crystal Ball function of m multiplied by a Gaussian func- 
tion of 5m. The Crystal Ball function models nearly all of the signal, including the 
radiative tail, while the Gaussian function models the effect of measurement errors in 
the direction transverse to m. In these fits, p, a, and as were allowed to float. The 
parameters p and a also appear in other terms; they float in all terms, but their values 
are constrained to be identical in all terms. The parameters a and n used in DT fits to 
data or Monte Carlo events with ISR, were fixed to the averages over all modes of the 
parameters obtained in ST fits to data or Monte Carlo events with ISR, respectively. In 
fits to Monte Carlo events without ISR the Crystal Ball function is replaced with a pure 
Gaussian. The values of Gaussian and Crystal Ball parameters obtained in DT yields are 
close to - but not exactly equal to - those obtained in the ST fits. 

• b(m; p, pl<7, Prv) b(m; fi, Pl&, Pro), the product of two bifurcated Gaussian functions. 
These functions account for misreconstructed signal events as they did individually in 
ST fits. The ratios pi and pR were obtained from the fits to ST Monte Carlo events 
without ISR. The signal contribution is the sum of this contribution and the Crystal Ball 
contribution. 

• a(m; Eq,£) g(5m; p, a c ), the product of an ARGUS background function of m and a Gaus- 
sian function of 5m. This term models continuum background and any residual mispar- 
titioned events remaining after choosing the value M(D) closest to M(D) in events with 
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multiple candidates. The Gaussian parameter a c is determined from continuum Monte 
Carlo events, and £ is allowed to float in the fit. 

• a(m; E , £ m ) g(fh] //, o^), the product of an ARGUS function of m and a Gaussian function 
of fa. This term models the horizontal band in the scatter plot. The values of £ m and 
<Jb are the identical in this term and the next. The value of o\, is determined from DT 
Monte Carlo events without ISR and is the same for all modes, while £ m is allowed to 
float in the fits. 

• a(fa; E Q , £ m ) g(m; fi, a^), the product of an ARGUS function of m and a Gaussian function 
of fa. This term models the vertical band in the scatter plot. The value of £ m in this 
term is constrained to be equal to the value of £ m in the term for the horizontal band. 

TABLE III: Double tag data yields and efficiencies and their statistical uncertainties. 



D Mode 




L> Mode 






Yield (10 2 ) 


Efficiency (%) 


D° - 


-> K-^ 




i) — > 




7T~ 




1.09 ±0.11 


42.6 ±0.5 


D° - 


-> K-^ 


-7T° 






7T~ 


TT° 


4.84 ±0.23 


12.1 ±0.3 


D° - 


-> K-^ 




D° — > 




7T~ 


7T~7r + 


2.80 ±0.17 


20.8 ± 0.4 


D° - 


-> K-^ 




5°^ 




7T~ 


^° 


2.45 ±0.16 


23.2 ± 0.4 


D° - 


-> K-^ 


-7T° 


5°^ 




7T~ 




2.62 ±0.16 


22.6 ±0.4 


D° - 


-> K-^ 




L> ^ 




7T~ 


7T~7r + 


2.05 ±0.14 


29.6 ± 0.4 


D° - 


-> K-^ 


- + — 
7T 7T 


5°^ 


K+ 


7T~ 




1.97 ±0.14 


29.6 ± 0.4 


D° - 


-> K-^ 






K+ 


7T~ 




3.59 ±0.20 


15.2 ±0.3 


D° - 


+ i^vr 4 




z)° — > 


K+ 


7T~ 


7T° 


3.40 ±0.19 


15.5 ±0.3 


D+ - 




f vr+ 




■K^ 


~7T~ 




3.79 ± 0.20 


26.7 ±0.4 


D+ - 






D" -s 


Kl 


.7T~ 




0.090 ± 0.030 


20.6 ± 0.4 


D+ - 






D~ -s 


Kl 


.7T" 




0.609 ± 0.079 


23.7 ±0.4 


D+ - 






D- -s 


■K^ 


"7T" 




0.530 ± 0.073 


23.9 ± 0.4 



Figures |3] and El show the projection of each of the "diagonal" M(D) vs. M(D) distributions 
projected onto the M(D) axes. The fit in each figure is the sum of all terms in the fit 
projected on the M(D) axis. Table ITTT1 gives the DT yields in data and the corresponding 
DT efficiencies calculated from the signal Monte Carlo events with ISR. These quantities 
are used in the fit - described in Sec. [X]for branching fractions and numbers of DD events. 

VIII. EFFICIENCIES IN DATA AND MONTE CARLO EVENTS 

We estimated efficiencies for reconstructing D decays using Monte Carlo simulations. For 
precision measurements, we must also understand the accuracy with which the Monte Carlo 
events simulate tracking efficiencies. Using partial reconstruction techniques, we measure 
tracking efficiencies - in both data and Monte Carlo events - for charged pions and kaons, 
as well as the efficiency for detecting K s and 7r° mesons. For example, we fully reconstruct 
one D° and then combine it with one or more of the other particles in the event to partially 
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reconstruct the D°, leaving out one particle - ir^, K ± , or K s - for which we wish to measure 
the efficiency. We calculate the square of the missing mass (MM 2 ) of the combination, which 
should peak at the mass squared of the missing particle. We then check whether the missing 
particle is actually detected or not, and we plot the missing mass squared for both classes 
of events. Fitting the signal peaks tells us what fraction of particles were actually detected. 
We also use similar techniques with %p' — > J/ip7c + ir~ and if)' — > <//■?/> 7r°7r° events to measure 
7r tracking efficiency at lower momenta and 7r° efficiency, respectively. 

We illustrate this technique for measuring K tracking efficiencies in Fig. |3 where - for 
both data and Monte Carlo events - we reconstruct D° mesons in the modes, D° — ► K~tt + , 
D° -> K-jt+it , and D° -> K-jt+it+tt-. We use the requirements \M(D°) - M D o\ < 
5 MeV/c 2 and \ AE\ < 25 MeV to obtain a clean sample of D° decays. Figures a) and^b) 
show the MM 2 distributions for data and Monte Carlo events, respectively, when the K ± 
was found, while Figs. Etc) and Efd) show the corresponding MM 2 distributions for events 
in which the K ± was not found. When the K ± is found, there is a clean peak with little 
background in both data and Monte Carlo events. Each of these distributions is fit with two 
Gaussians of the same mean but different widths, and a small constant background. The 
widths of the Gaussians are somewhat different in data and Monte Carlo events. The MM 2 
distributions for events in which the K was not found are more complex. There is a peak 
at M 2 + due to D° — > decays, and a peak at M^+ due to D° — > K + n~ decays. The 

shoulder on the right of each figure is due to D° — > A + 7r~7r° decays in which only the ir~ 
is detected; its shape is well approximated by an error function. Finally, there is a linearly 
rising shape which describes the remaining backgrounds. We determine the parameters for 
the shapes of these backgrounds in fits to signal Monte Carlo events for these processes. We 
fit the data in Figs. Efc) and E^d) using these background shape parameters and the widths 
of the Gaussian parameters determined in fits to events in which the was detected. 

For ^ and K^, we find evidence for slightly higher tracking efficiencies in data than 
in Monte Carlo events. Based on these preliminary studies, we increase the efficiencies 
determined from the Monte Carlo simulations by a factor of 1.03 ± 0.03 for each charged 
track candidate used in the event. In our preliminary studies of 7r° efficiencies, we see no 
evidence for an efficiency difference between data and Monte Carlo events. Also, we find 
no evidence for a correction for Kg decays beyond the factor of 1.03 2 for the two charged 
daughter pions in the decay. 

IX. SYSTEMATIC ERRORS 

We take systematic uncertainties into account directly in the branching fraction fit. Ta- 
ble |TV] lists the uncertainties that we included and a brief description of each contribution 
follows. 

• Data processing 

Potential losses of events in data processing has been limited to 0.3% and is assessed by 
reanalyzing the data with slightly different software configurations. 

• Yield fit functions 

We gauge the sensitivity of the ST and DT yields to variations in the M(D) fit functions 
by repeating the fits with different values of pa, Pr, and pi as well as with different 
functional forms (such as a symmetric Gaussian instead of a bifurcated Gaussian for the 
wide signal component) and comparing efficiency-corrected yields. 
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TABLE IV: Systematic uncertainties and the quantities to which they are applied in the branching 
fraction fit. Uncorrelated, additive uncertainties are given in the first section, and correlated, 
multiplicative uncertainties in the second. 



UUUl i^C 


htp r*i"i fin 5i 1 T In r'P'T'l" m - n1"\7' ( % i 






u.o 


All tritilrlG 


1 ItlU. lllj 1 U.11U LlUllo 


n 1-9 q 


All viplH<3 

-riii yiciLio 


ri o pLrrrrrii inn Kiqc 
JDdL-K^l UUI1U. Uldo 


z .0 


LJ ± y lclU-o 


L/UUUlc LJKjOU lIlLcl lei cllCc 


u.o 


TVfin+ral FIT viplrlQ 
IV cULl dl U A. yiclLlo 


iyCUCOlul ollllU-lcLtlCIl 


o . u 


TrPi pl/i ti cr tifii piphpipq 




3.0 


Ko efficiencies 




4.4 


7T° efficiencies 




0.3 


■k^ PID efficiencies 




1.0 


PID efficiencies 


Trigger simulation 


0.3 


ST efficiencies 


Final state radiation 


0.5 


D efficiencies 


\AE\ requirement 


1.0 


D efficiencies, correlated by decay 


Resonant substructure 


3.0 


D° — > i^ _ 7T + 7r + 7r~ efficiencies 



• Background bias 

We estimate possible biases in the background parametrization from the agreement of the 
fit functions with the data. We find the discrepancies to be negligible for ST yields and as 
large as 2.5% for DT yields. Additional bias in the DT yields could have been introduced 
by the procedure for selecting the best candidate per event in double tags. However, as 
described in Sec. IV111 we searched for this type of background in simulated continuum 
and generic DD Monte Carlo events and found no evidence of such a contamination. 

• Double DCSD interference 

In the neutral DT modes, the Cabibbo-allowed amplitudes can interfere with amplitudes 
where both D° and D° undergo doubly Cabibbo-suppressed decays (DCSD). We include 
an uncertainty to account for the unknown phase of this interference. 

• Detector simulation - tracking, K s , and ir° efficiencies 

We estimate uncertainties due to differences between efficiencies in data and those esti- 
mated in Monte Carlo simulations using the partial reconstruction technique described 
in the previous section. 

• Detector simulation - particle identification (PID) efficiencies 

Particle identification efficiencies are studied by reconstructing decays with unambiguous 
particle content, such as D° — > K s 7t + 7t~ and <fi — > K + K~ . We also use D° — > K~tt + 7t , 
where the K~ and tt + are distinguished kinematically. The efficiencies in data are well- 
simulated by the Monte Carlo, and we assign correlated uncertainties of 0.3% and 1.0% 
to each n + and K + , respectively. We do not assign these uncertainties to K s daughters, 
because they are not subjected to the n + identification requirements. 
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• Trigger simulation 

We assign a conservative uncertainty in the trigger simulation efficiency taken from the 
trigger inefficiency we find in simulations of the D° — > K~-k + tt mode. 

• Final state radiation 

In Monte Carlo simulations, final state radiation typically reduces DT efficiencies by a 
factor of approximately two times the reduction of ST efficiencies. This leads to branching 
fraction values larger by 0.5% to 2% than they would be without including FSR in the 
Monte Carlo simulations. We have verified the accuracy of the FSR simulation to roughly 
10% of itself using ip(2S) — ► J/tfmn decays, where J/ip — > fi + fi~-y. For these preliminary 
results, we assign conservative uncertainties of 0.5% to ST efficiencies and 1.0% to DT 
efficiencies, correlated across all modes. This results in a ±0.5% correlated systematic 
error in branching fractions. 

• | AS | requirement 

Discrepancies in detector resolution between data and Monte Carlo simulations can result 
in differences between the efficiencies of the AE requirement in data and Monte Carlo 
events. No evidence for such discrepancies has been found, and we include a systematic 
uncertainty of 1.0% for ST modes and 2.0% for DT modes, correlated by decay process. 

• Resonant substructure 

The resonant substructure of D° — > K~7t + 7t + 7t~ in data is found to disagree with that 
in Monte Carlo simulations, introducing a possible fractional bias in efficiency of 3.0%. 
We apply a correlated systematic uncertainty of this size to all ST and DT modes with 
a D° — ► K~7t + 7t + 7t~ decay. 

For the e + e~ — ► DD cross section measurements, we include additional uncertainties from 
the luminosity measurement (3.0%) and from variations in the center-of-mass energy. We 
estimate the latter uncertainty by calculating N D oi)0 and N D + D - in data taken at different 
values of E cm in a 4 MeV region on the ^(3770). We find fractional variations of 2.5% in 
N d oqo and 1.4% in N D + D -, and include these numbers as systematic uncertainties. 

Finally, we note that - although there are many contributions to the systematic errors 
listed in Table IIVI - the contribution of tracking efficiencies dominates in this preliminary 
result. This is due to our conservative assignment of an error of ±0.03n, where n is the 
number of tracks in a ST or DT mode for this uncertainty. 



X. RESULTS OF THE BRANCHING FRACTION FITS 

To determine the branching fractions B(D° -> IC'ir*), B(D° -> R-tv+tt ), B(D° -> 
K-n + n + n-), B{D + -> K-n + n + ), and B(D + -> K° s n + ) as well as N D o D o and N D+D -, 
we measure event yields and efficiencies for the ten ST modes and thirteen DT modes, 
given in Tables ITD and UTTl In the branching fraction fit, we correct these event yields not 
only for efficiency but also for crossfeed among the ST and DT modes and backgrounds 
from other D decays. The estimated crossfeed and background contributions induce yield 
adjustments of 2% at most. Their dependence on the fit parameters is taken into account 
both in the yield subtraction and in the x 2 minimization. In addition to the correlated and 
uncorrelated systematic uncertainties, the statistical uncertainties on the yields, efficiencies, 
and background branching fractions are also included in the fit. 

We validated the algorithm and the performance of the branching fit - as well as our 
entire analysis procedure - by measuring the branching fractions in generic Monte Carlo 
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events. Since the generic Monte Carlo simulations were done without ISR, we used effi- 
ciencies obtained from signal Monte Carlo events without ISR. We find that the results of 
this procedure were in excellent agreement with the branching fractions used by EvtGen in 
generating the events. Furthermore, the generic Monte Carlo sample has an order of mag- 
nitude more events than our data, so the statistical errors in this test were about a factor 
of three smaller than the errors in the branching fractions obtained from data. The system- 
atic errors were also substantially smaller than those estimated for data, so the agreement 
between measured and generated branching fractions of the generic Monte Carlo events is a 
very stringent test of our entire analysis procedure. 



TABLE V: Preliminary fitted branching fractions and DD pair yields. Uncertainties are statistical 
and systematic, respectively. 



Parameter 




Fitted Value 


Fractional 


Fractional 










Stat. Error 


Syst. Error 


N D o D o 






(1.98 ±0.04 ±0.03) x 10 5 


1.9% 


1.5% 


B(D° - 


> 




0.0392 ± 0.0008 ± 0.0023 


2.1% 


5.8% 


B(D° - 




-7T ) 


0.143 ±0.003 ±0.010 


2.0% 


7.2% 


B(D° - 






0.081 ± 0.002 ± 0.009 


2.1% 


11.6% 


N D+D - 






(1.48 ± 0.06 ± 0.04) x 10 5 


4.3% 


2.4% 


B(D + - 


-> K~ir' 




0.098 ± 0.004 ± 0.008 


4.4% 


8.5% 


B{D+ - 




-) 


0.0161 ±0.0008 ±0.0015 


4.9% 


9.1% 


B(D° - 




-7T°)/B(D -> K~ir+) 


3.64 ±0.05 ±0.17 


1.3% 


4.8% 


B(D° - 




-7T + Tr-)/B{D° -» K-7T+) 


2.05 ±0.03 ±0.14 


1.4% 


6.8% 


B(D + - 




~)/B(D + -» i^-7r+vr+) 


0.164 ±0.004 ±0.006 


2.5% 


3.7% 



The results of the data fit are shown in Table IVl The \ 2 °f the fit is 8.9 for 16 degrees of 
freedom, corresponding to a confidence level of 92%. To obtain the separate contributions 
from statistical and systematic uncertainties, we repeat the fit without any systematic inputs 
and take the quadrature difference of uncertainties. All five branching fractions are consistent 
with - but also higher than - the current PDG averages If no FSR is included in the 
simulations to calculate signal efficiencies, then all the branching fractions would be 0.5% to 
2% lower. In addition, we compute ratios of branching fractions to the reference branching 
fractions, which have higher precision than the constituent branching fractions, and these 
also agree with the PDG averages. These ratios are also sensitive to final state radiation, 
and - without these corrections - all three ratios would be 1% to 2% higher. 

The correlation matrix for the seven fit parameters is given in Table IVII In the absence 
of systematic uncertainties, there would be no correlation between the charged and neutral 
D parameters. That these correlations are, in fact, large indicates that the precision of our 
branching fraction measurements is limited by our understanding of systematic effects. 

We obtain the e + e~ — > DD cross sections by dividing the fitted values of N d odo and 
Nd+d- by the luminosity collected on the ?/>(3770), which we determine to be C = 57.2 ± 
1.7 pb -1 . Thus, at E cm = 3773 MeV, we find preliminary values of the production cross 
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TABLE VI: Correlation matrix, including systematic uncertainties, for the branching fractions 
determined from the fit. 





at - i*r~Tr+ tr~ -T7-+-T7-0 tf- 

JV£)0£)0 J\ TT A l 1 A 




AT , 

^d+d- 




K S TT 


-/V no nO 


1 -0.37 -0.26 
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-0.02 


-0.02 


K--K+ 


1 0.75 


0.90 


0.00 


0.76 


0.69 


K~TT + TT° 


1 


0.73 


0.00 


0.62 


0.57 


K~ir + TT~1T + 




1 


0.00 


0.80 


0.73 


N D+D - 






1 


-0.49 


-0.47 


K~TT + TT + 








1 


0.90 
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sections, 

a(e + e- -> D°D°) = (3.47 ± 0.07 ± 0.15) nb 

a(e + e- -> D+D~) = (2.59 ± 0.11 ± 0.11) nb 

a(e + e~^DD) = (6.06 ± 0.13 ± 0.22) nb 

a(e + e- ^ D + D-)/a(e + e- ^ D°D°) = 0.75 ± 0.04 ± 0.02, 

where the uncertainties are statistical and systematic, respectively. The charged and neutral 
cross sections have a correlation coefficient of 0.32 stemming from the common use of the 
luminosity measurement. 



XI. CONCLUSIONS 

We report preliminary measurements of three D° and two D + branching fractions in a 
sample of nearly 60 pb _1 of e + e~ — > DD data obtained at E cm = 3.77 GeV. The results are 
presented in Table IVl and the correlation coefficients among the results are given in Table IVTl 
We find branching fractions in agreement with - but somewhat higher - than those in the 
PDG 0] summary. We note that the Monte Carlo simulations used in calculating efficiencies 
in this analysis included final state radiation. However, the branching fractions used in the 
Particle Data Group averages do not include this effect. If we had not included final state 
radiation in our simulations, branching fractions would have been 0.5% to 2% lower. At this 
preliminary stage in our analysis, the systematic errors are factors of approximately 2 to 5 
larger than the statistical errors. Our statistical errors for branching fractions and ratios of 
branching fractions are - in nearly all cases - less than the statistical errors in the individual 
measurements that contributed to the PDG averages. 

Our techniques for event selection, for fitting data and Monte Carlo distributions to 
obtain yields and efficiencies, and for fitting the results to obtain branching fractions are 
validated with generic Monte Carlo events to a much higher precision than the statistical er- 
rors from the current data sample. This indicates that we can improve the results from these 
data substantially with more detailed understanding of any possible systematic differences 
between our data and our Monte Carlo simulations. Our systematic errors are dominated by 
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the uncertainty of 0.03n, where n is the number of charged tracks in a ST or DT mode, that 
we assign for the difference between tracking efficiency in data and Monte Carlo simulations. 
Many of the systematic uncertainties, such as those for tracking and particle identification 
efficiencies, will be improved with larger data samples. 

Our preliminary measurements of the production cross sections a(D°D°), a(D + D~), 
and a(DD) are in good agreement with the cross sections measured by the MARK III 
collaboration and in excellent agreement with recent BES results reported at the 2004 
Moriond conference (l7| . 
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FIG. 1: Distributions of calculated M(D) values for single tag neutral D candidates in data with 
D° and D° combined in each mode. The points are data and the curves are fits to the data. 
In each plot, the dotted curve shows the contribution of the ARGUS background function, the 
dashed curve shows the contribution of the bifurcated Gaussian, and the solid curve shows the 
total, including the Crystal Ball function, (a) and (b) illustrate D° — > K~ir + candidates, (c) and 
(d) illustrate D° — > K~ir + iT candidates, and (e) and (f) illustrate D° — > K~ir + 7r + ir~ candidates. 
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FIG. 2: Distributions of calculated M(D) values for single tag charged D candidates in data with 
D + and D~ combined in each mode. The points are data and the curves are fits to the data. 
In each plot, the dotted curve shows the contribution of the ARGUS background function, the 
dashed curve shows the contribution of the bifurcated Gaussian, and the solid curve shows the 
total, including the Crystal Ball function, (a) and (b) illustrate D + — > K~ir + ir + candidates and 
(c) and (d) illustrate D + — > K$ir + candidates. 
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FIG. 3: Scatter plot of M(D) vs. M(D) values for D° -» K+tt-tt vs. D° -► R-it+tt DT 
candidates from simulated events, (a) illustrates all candidates after imposing the AE requirement 
for each candidate given in Table QJ (b) illustrates the candidates remaining after choosing the 
candidate with the average M{D) closest to M D o in events with multiple candidates. 
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FIG. 4: Projections of double tag D°-D° candidate masses on the M(D) axis, with the D° and 
D° reconstructed in corresponding charge conjugate modes. Projections of the candidate masses 
on the M(D) axis are nearly identical, (a) illustrates D° — > K~tt + candidates, (b) illustrates 
D° — > K~ir + ir° candidates, and (c) illustrates L>° — » i"T~7r + 7r + 7r~ candidates. 
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FIG. 5: Projections of double tag D + -D~ candidate masses on the M(D) axis, with the D + 
and D~ reconstructed in corresponding charge conjugate modes. Projections of the candidate 
masses on the M(D) axis are nearly identical, (a) illustrates D + — > K~ir + 7r + candidates and (b) 
illustrates D + — > K%n + candidates. 
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FIG. 6: Missing mass squared distributions for measuring efficiencies in data and Monte Carlo 
events in which a D° has been found, (a) and (b) data and Monte Carlo events, respectively, in 
which the K + from D° — > K + ir~ decay was found, (c) and (d) data and Monte Carlo events, 
respectively, in which the K + from D° — > K + n~ decay was not found. 
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